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i The diffractive production of vector mesons in deep inelastic scattering (DIS) is calcu- 

. lated in the McLerran-Venugopalan model. This is relevant when large parton densities 

are probed by the virtual photon, as is the case at small Bjorken x or in DIS off nuclei. 
We investigate differences between the exclusive production (when the target doesn't 
break up) which dominates at small momentum transfer squared \t\, and the diffractive 
production (when the target scatters inelastically) which dominates at large 

£h! 1 Motivation 

' Diffractive vector meson production in deep inelastic scattering (DIS) ~f*A — > VY, where A 

i-C ■ stands for the target nucleus and Y for the final state it may dissociate into, is a process 

in which elastic and inelastic interactions of the target can be experimentally distinguished. 
t-H ■ At high energies, the qq dipole that the virtual photon has fluctuated into scatters off the 

gluonic field of the nucleus before recombining into the vector meson. While this scattering 
involves a color-singlet exchange, leaving a rapidity gap in the final state, the nucleus can 
still interact elastically (Y = A, this is called coherent diffraction) or inelastically (i.e. break 
t^j- ■ up, called incoherent diffraction). 

| Kinematically, a low invariant mass of the system Y corresponds to a large rapidity gap in 

the final state between that system and the vector meson, and implies that the longitudinal 
' momentum of the meson is close to that of the incoming photon. In this case, the eikonal 

approximation can be assumed to compute the dipole-nucleus scattering. At small values of 
x = (Q 2 + My)/(Q 2 + W 2 ) where Q 2 is the photon virtuality, My the vector meson mass, 
• *h . and W the energy of the 7* — A collision, a target proton can also be considered. Indeed in 

that case, since partons with an energy fraction as small as x are probed in the target wave 
function, the dipole will scatter off large gluon densities generated by the QCD evolution. 

The cross-section is maximal at minimum momentum transfer with exclusive production 
(or coherent diffraction) dominating. As the transfer of momentum gets larger, the role of 
incoherent diffraction increases and eventually it becomes dominant, typically for momenta 
bigger that the inverse target size; the elastic contribution decreases exponentially while the 
inelastic contribution decreases only as a power law. In the case of a target proton, it is 
known that saturation models describe well the exclusive cross section [2J, while the BFKL 
Pomeron exchange approach works well for the target-dissociation cross-section [3] . 

In this work we show that, with the Color Glass Condensate (CGC) picture of the 
small— 2; part of the hadronic/nuclear wave function, both coherent and incoherent diffrac- 
tion can be described in the same framework, for protons and nuclei. We also explicitly 
calculate both contributions to the diffractive vector meson production cross-section using 
the McLerran-Venugopalan (MV) model for the CGC wave function. Finally, we discuss 
phcnomenological consequences in the context of a future electron- ion collider [4j. 
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2 Diffractive Vector Meson Production in DIS off the CGC 



The CGC is an effective theory of QCD [S] which aims at describing the small x part of the 
hadronic/nuclear wave function, when the gluon density is so large that non- linear effects 
are important. Rather than using a standard Fock-state decomposition, it is more efficient 
to describe it with collective degrees of freedom, more adapted to account for the collective 
behavior of the small- x gluons, which do not interact with a probe independently, but rather 
behave coherently. The CGC approach uses classical color fields: 

\h) = \qqq) + \qqqg) +••• + \qqqg ■ ■ ■ ggg) + ■■■ => \h) = J DA <S> XA [A] \A) . (1) 

The long-lived, large- a; partons are represented by a strong color source p which is static 
during the lifetime of the short-lived small- x gluons, whose dynamics is described by the 
color field A^ l/.9s- X A denotes an arbitrary separation between the field and the source. 

The CGC wavefunction <& XA [A] is the fundamental object of this picture, it is mainly a 
non-perturbative quantity, but the xa evolution can be computed perturbatively. Requiring 
that observables are independent of the choice of xa, a functional rcnormalization group 
equation can be derived. In the leading-logarithmic approximation which rcsums powers of 
as \n(l/xA) 7 the JIMWLK equation [6] describes the evolution of [^4] | 2 with xa- The 
information contained in the wavefunction, on gluon number and gluon correlations, can be 
expressed in terms of n-point correlators, probed in scattering processes. These correlators 
consist of Wilson lines averaged with the CGC wavefunction, and resum powers of gsA. 

In diffractive vector meson production, the relevant quantity is (the photon is a right 
mover, the CGC a left mover, and the gauge is A + = 0): 



T xy [A-} = 1 - ^-Tr (U y U x ) , with U X [A~] = V exp (ig s J dz+T c A- (z+ 



(2) 



In terms of this object, the differential cross sections for a transversely (T) or longitudinally 
(L) polarized photon are given by (with t = — q\ the momentum transfer squared) 



dt 



J dzd 2 xd 2 ye lq ^ ■ ("K 1 "*)^ * T>L (z, x- y)T xy 



(3) 



where 2* T = + ^>~~ and *l = with 



^(z,r) = ^[#(z,r)]>f (z,r) , (4) 



hh 



the overlap between the photon and meson wave functions. A and h denote polarizations 
and helicitics while z is the longitudinal momentum fraction of the photon carried by the 
quark and x and y are the quark and antiquark positions in the transverse plane. 
The target average ( . ) x is done with the CGC wave function squared l^^^ - ] | 2 : 

(/)» = / DA-\$ x [A-\\ 2 f[A-\ . (5) 

If one had imposed elastic scattering on the target side to describe the exclusive process 
-f*A — > VA, the CGC average would be at the level of the amplitude, and the two-point 
function (T xy ) x inside the | . | 2 in ([3]), recovering the formula often used with dipolc models. 
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3 The McLerran-Venugopalan model 

Instead, when also including the target-dissociative part, the diffractive cross section involves 
the 4-point correlator (T xy T uv ) x . In order to compute it, we must specify more about the 
CGC wave function. We shall use the McLerran-Venugopalan (MV) model [7], which is a 
Gaussian distribution for the color charges which generate the field A : 

l*,L4-]| 2 = cxp (-J #x<Pydz+ Pc{z+ >*^* + ' y) ^ , (6) 

where the color charge p c and the field A~ obey the Yang-Mills equation — V 2 *4~ (z + , x) ~ 
gsp c (z + ,x). The variance of the distribution is the transverse color charge density squared 
along the projectile's path p 2 (z + ), with 

( Pc (z + ,x) Pd (z' + ,y)) = 6 cd 5(z+~z'+)6^( X -y)p 2 (z+) . (7) 

The only parameter is the saturation momentum Q s , with Q 2 proportional to the integrated 
color density squared. Note that there is no x dependence in the MV model, it should be 
considered as an initial condition to the small— x evolution. 

The MV distribution is a Gaussian distribution, therefore one can compute any target 
average by expanding the Wilson lines in powers of gsA~ (see @), and then use Wick's 
theorem [5] . All correlators of A's can be written in terms of 

g%{A- (x+, x)A d (y+,y)) = S cd 5(x+ - y+)g 4 s p 2 (x + ) J d 2 z G(x - z)G(y - z) , (8) 
with the two-dimensional massless propagator 

|k|>A QC D 

The color algebra is the difficult part to deal with. 

We give the result for the (S xy S uv } correlator [9], with S xy = 1 — T xy : 



/ Q q \_ „-^f [F(x-y)+F(u-v)] 

= (S X y)(S uv ) 

( F(x,u;y,v)-y/A _ F(x,y; u,v) \ 
\ 2^A N 2 VA 



F(x, u;y,v) + VA _ F(x,y;u,v) ^ ^^^A 

2^A n?Va 



e 4 



(10) 

where p 2 = J dz + p 2 (z + ). This seemingly complicated object is given in terms of a single 
F(x - y) = gtp 2 J d 2 z [G(x - z) - G(y - z)] 2 . (11) 



function: 



Indeed also one has 



A = F 2 (x, u; y, v) + F(x, y; u, v)F(x, v; u, y) , (12) 

^ C 

-2 M 2 F(x,y;u,v) = F(x-u)+F(y-v)-F(x-v)-F(y-u) . (13) 
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In the function F(r), the infrared cutoff Aqcd only enters through a logarithm as ex- 
pected. In the |r|Agc , £)^;l limit, one has 



2 W 2tt h 



dk 



QCD 



1 - Jo(fc|r|) 

fc 3 



r 

T 



An 



H 2 log 



r 2 A2 

1 QCD 



(14) 



This is the standard definition of the saturation scale in the MV model. It is also possible 
to consistently include small— a; evolution in a such a calculation, in the large— N c limit [lOj 
and beyond [11]. Essentially one should replace n(z + ) by fi x (z + ,x,y) in ([6]), and the x 
evolution of the corresponding F(x, x, y) can be obtained from the JIMWLK equation. 

Finally, note that in (fT0|) for x = y or u = v, F(x, y; u, v) = and one recovers the single 
dipolc average. Also, in the large— N c limit, one has (T X yT uv ) = (T xy )(T uv ), which means 
that at small— x, the target-dissociative part of the diffractive cross-section in suppressed at 
large 2V C , compared to the exclusive part. 



4 Results 



The results presented in this section are obtained 
with the x evolution of the saturation scale mod- 
eled as in [T2l: 



■w - (?) 



A/2 



GeV , 



(15) 



with A = 0.277 and x = 4.1 10~ 5 for the case of 
a target proton. The collincar logarithm of Q s 
(see (O) is neglected, which corresponds to ex- 
act geometric scaling [13]: F(x,r) = F[r 2 Q 2 (x)]. 
As an illustration, the resulting cross-section for 
diffractive J/^> production is displayed in Fig.l, 
and separated into its coherent and incoher- 
ent contributions. The light-cone Gaussian J/'I' 
wave function [14] has been used in (j4j. At small 
values of \t\ where coherent diffraction domi- 
nates, our results are in agreement with HERA 
data [TS] (one can get a better agreement with 
more realistic saturation models [2], but this is 
not our point). Our model indicates that for 
\t\ > 0.7 GeV 2 or so (this value slightly decreases 
when Q 2 increases), incoherent diffraction starts 
to dominate. This may be the reason why the 
data on exclusive production stop: there is too 
much proton-dissociative 'background'. We ob- 
serve that this part of the cross-section decreases 
as a power law with \t\, rather than exponentially 
as the exclusive part does. 




Figure 1: Diffractive J/^P production in 
DIS at HERA, for W = 90 GeV and dif- 
ferent values of Q 2 . In our calculation, we 
separated the exclusive (full lines) and the 
proton-dissociative (dashed lines) parts. 
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In the case of a target nucleus, we expect the following qualitative changes in the t 
dependence. First, the low— \t\ regime with elastic scattering of the nucleus will be dominant 
up to a smaller value of \t\ compared to the proton case, reflecting the bigger size of the 
nucleus. Then, the nucleus-dissociative part will be split into two: an intermediate regime 
in momentum transfer up to about 0.7 GeV 2 where the nucleus will predominantly break up 
into its constituents nucleons, and a large — 1£| regime where the nucleons inside the nucleus 
will also break up, implying pion production in the Y system for instance. 

The model discussed in this work is well adapted to describe the low- and large — 1£| 
regimes, but not the intermediate regime since the constituent nucleons are absent from 
the description JTJ). This problem has been addressed in a complementary setup in the 
case of inclusive diffraction off nuclei |16j . and the coherent diffraction regime was found 
to be dominant up to about \t\ = 0.05 GeV 2 . The vector meson production case will be 
addressed next. While in the proton case, both exclusive and diffractive processes can be 
measured, it is likely that at a future electron-ion collider, the exclusive cross section cannot 
be extracted: when the momentum transfer is small enough for the nucleus to stay intact, 
then it will escape too close to the beam to be detectable. Therefore the diffractive physics 
program will rely on our understanding of incoherent diffraction. 
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